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We present an experimental and theoretical analysis of the self- and cross-Kerr effect of extended 
plasma resonances in Josephson junction chains. We calculate the Kerr coefficients by deriving 
and diagonalizing the Hamiltonian of a linear circuit model for the chain and then adding the 
Josephson non-linearity as a perturbation. The calculated Kerr-coefficients are compared with the 
measurement data of a chain of 200 junctions. The Kerr effect manifests itself as a frequency shift 
that depends linearly on the number of photons in a resonant mode. By changing the input power 
on a low signal level, we are able to measure this shift. The photon number is calibrated using the 
self-Kerr shift calculated from the sample parameters. We then compare the measured cross-Kerr 
shift with the theoretical prediction, using the calibrated photon number. 

BAGS numbers: 74.50.+r,74.81.Fa,85.25.Cp,42.65.Hw 


I. INTRODUCTION 

One-dimensional arrays of Josephson junctions, or 
Josephson junction chains, have received considerable in¬ 
terest for more than three decades. They were originally 
introduced as a theoretical model system for the study 
of the zero-temperature superconductor-insulator tran¬ 
sition in superconducting granular filmsii^. Josephson 
junction chains consist of superconducting islands with a 
small capacitance to ground Co, connected to each other 
by Josephson junctions. The superconductor-insulator 
transition originates from the Josephson potential energy 
Ej cos(0) as this allows the winding by 2tt of the phase 
difference (f) between neighboring islands. At zero tem¬ 
perature, these phase-slips are driven by the quantum 
fluctuations of the phase induced by Coulomb charging 
effects. As a result, the transition occurs as a function of 
the parameter EjjEo, where Eq = I2Cq the Coulomb 
charging energy of the island. Since the quantum fluctu¬ 
ations of the phase actually correspond to the zero-point 
fluctuations of propagating electromagnetic modes along 
the chain^, the transition can also be characterized in 
terms of the impedance Zj Rq = ^E^^jEj associated 
with the propagation of these modes, where Rq is the 
superconducting resistance quantum Rq = hjAe^. In¬ 
deed, the transition to the insulating state occurs when 
ZjRQ becomes larger than ^ 7r/2. The superconductor- 
insulator transition has been observed in granular films^ 
and wires^, as well as in Josephson junction chains^i^. 

More recent work focused on Josephson junction chains 
of finite length for which the junction capacitance C is 
much larger than the capacitance Cq to ground. In such 
chains, the above mentioned 27r phase windings occur in 
the form of coherent quantum phase-slips (QPSs). It has 
been shown^“— that coherent QPSs occurring locally at 
each junction of the chain induce a nonlinearity related 
to the global charge q on the chain Uo cos(7rg/e), dual to 
the usual Josephson nonlinearity. Here, the amplitude 


Uo = Nuo scales with the length N of the chain. The 
amplitude uq {E^EcY^^ exp — yj 8Ej / Ec is the ampli¬ 
tude for a single QPS to occur on one of the junctions of 
the chain, where Ec = e^/2C. In this limit, the chain be¬ 
haves as a so-called quantum phase-slip junction (QPSJ), 
a device dual to the usual Josephson junction and whose 
properties have been discussed since the pioneering work 
by Likharev and coworker o^^d^ . For instance, a QPSJ 
is expected to sustain Bloch oscillations, dual to Joseph¬ 
son oscillations. Evidence of Bloch oscillations has been 
found in Refs. [IMi The actual realization of a QPSJ 
could have far-reaching consequences in quantum metrol¬ 
ogy and quantum information processingi^^— . 

The effect of the propagating modes of the chain on the 
properties of the QPSJ has not been taken into account 
in Refs. This is correct for relatively short chains, 

with N < 2C/Co. Yet for longer chains, these modes 
appear and their effect on the phase-slip amplitude is 
non-negligibl o^^dQ . Indeed, in the thermodynamic limit, 
it is the interplay between the modes and the phase-slips 
that leads to the superconductor-insulator transition. In 
view of possible QPSJ applications, it is therefore impor¬ 
tant to study the modes directly, in the presence of the 
chain’s nonlinearity. This is the subject of the present 
paper, where we present two-tone spectroscopy measure¬ 
ments of the modes for a chain containing 200 junctions 
and tuned in the weakly nonlinear regime. We show that 
the results can be interpreted in terms of a model Hamil¬ 
tonian that takes the weak nonlinearity into account via 
self-Kerr and cross-Kerr interaction terms of the propa¬ 
gating modes. 

The paper is organized as follows. In Section [III we 
introduce the sample and the experimental setup and 
show spectroscopy measurements of the harmonic modes 
of a chain of 200 junctions. In Section illll we derive 
the Hamiltonian of a Josephson junction chain and its 
eigenvalues and eigenvectors in the linear limit. We also 
introduce the lowest order nonlinear terms as a pertur- 
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bation to the linear Hamiltonian. The two most promi¬ 
nent terms are the self-Kerr and cross-Kerr terms that 
appear. In Section lEl we first use the measured self- 
Kerr frequency shift to calibrate the photon numbers in 
the resonant modes of the Josephson junction chain. We 
then use this calibration to measure the cross-Kerr shifts. 
We find cross-Kerr shifts of the same order as the ones 
predicted by theory. In the last section, we summarize 
our results and conclude. 


II. EXPERIMENT 


A. Sample 

We have measured the resonances associated with 
the propagating modes of a Josephson junction chain 
containing 200 SQUIDs (Superconducting Quantum 
Interference Device). The sample is cooled to 20mK 
in a dilution refrigerator. The resonant modes are mea¬ 
sured with a vector network analyzer (VNA). A scheme 
of the measurement setup is drawn in figure [TJi. The 
input line is attenuated by —20dB at 4K and —AOdB 
at base temperature. The input bandwidth is limited by 
the coax-cables which are specified up to 18 GHz but still 
show transmission up to higher frequencies. The output 
line of the measurement setup comprises two cryogenic 
broadband isolators (Pamtech CW1019-K414), a cryo¬ 
genic amplifier (Caltech Cryo 1-12 SN262D) and three 
room temperature amplifiers. The output line has a gain 
of about 85 dB over a frequency band between 4 GHz and 
12 GHz. 

A SEM-image of the sample is shown in figure [TJd. The 
chain is coupled to 50 Cl microstrip transmission lines at 
both ends through coupling capacitors. In figure [it we 
show a zoom on the input coupling capacitor. Input 
and output capacitor have different sizes, so the chain 
is coupled to the output line stronger than to the input 
line. Junctions and transmission lines are fabricated in 
the same fabrication step by shadow evaporation of alu¬ 
minum on a silicon substrate. We use 100 keV electron 
beam lithography with an asymmetric undercut in the 
two-layer resist (PMMA/MAA, PMMA) to deposit all 
unwanted structures on the resist wall o^^i^^ . 

The sample parameters are summarized in Table HI The 
Josephson energy was determined by measuring the re¬ 
sistance of test junctions on the same chip at room tem¬ 
perature and by assuming an increase of the resistance of 
30% during cooldown, which is what is typically observed 
for our samples. The junction capacitance is determined 
from the junction size obtained by SEM-observation of 
the test junctions. We use a capacitance to surface ratio 
of 48fE//im2^i24. 



FIG. 1. (a) Schematic of the measurement setup, (b) SEM- 
image of the sample. On the top and bottom side one can 
see the ends of the microstrip transmission lines. The SQUID 
chain in between is barely visible on this scale, (c) Zoom on 
one of the coupling capacitors. The image is rotated 90° with 
respect to the top image. The insets show zoom on the chain 
and on a single SQUID. 


Ej (GHz) 

Ec (GHz) 

^ (GHz) 

G (fF) 

Lj (nH) 

40.1 ±2 

4.7 ±0.5 

38.4 ± 3 

4.2 ±0.4 

4.1 ±0.2 

Ic (uA) 

C„(fF) 

C'ont(fF) 

CICo 

Co (aF) 

81 ±4 

8±2 

24 ±6 

35 ±3 

120 ± 15 


TABLE I. Device parameters of the sample. 


B. Measurements 

Three of the modes of the chain lie within the band¬ 
width of the amplification chain and can be observed di¬ 
rectly by applying a single microwave tone and measuring 
the transmission coefficient S 21 between ports 1 and 2 of 
the VNA. The frequencies of the other modes lie outside 
the measurement bandwidth. 

As an example for the direct spectroscopy of reso¬ 
nances within the experimental accessible bandwidth, we 
plot the transmission magnitude 5'2i in figure [2] for the 
mode corresponding to j = 3. The spectrum is divided 
by a reference spectrum recorded at half flux frustration 
of the SQUID loop. In this way, the transmission back¬ 
ground is normalized to 1. Due to a parasitic coupling 
between the two transmission lines (input and output 
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FIG. 2. Transmitted amplitude through the j=3 resonance. 
The red curve shows a fit with the Fano formula eq. m 
using the following fit parameters: yo = 0.89, q = —3.17, 
j — 77.3 MHz and /o = 7.91 GHz. The inset shows the trans¬ 
mitted phase as a function of frequency for the same reso¬ 
nance. 



FIG. 3. Two-tone spectroscopy: phase of the transmission 
through the j = 2 resonance as a function of the frequency 
of an additional probe tone applied by an external microwave 
source. We observe the lowest 14 resonant modes of the chain. 
The blue datapoints at higher frequencies were measured with 
higher excitation power of the probe tone. 


lines) the resonances show a Fano line shap o^^i^^i^^ . We 
thus fit the resonances with the Fano formulai^ for the 
transmitted amplitude ^(/), 


y{f) = yo 


(g + 2(/-/o)/7)^ 

1 + 4(/-/o)V72‘ 


( 1 ) 


Here /o is the resonance frequency, 7 the width of the 
resonance (full width half maximum FWHM), yo is the 
amplitude related with the parasitic coupling between 
the two transmission lines, and q is the Fano factor that 
is given by the ratio between the amplitudes transmitted 
at resonance through the chain and the parasitic trans¬ 
mission. 

Modes outside the band can nevertheless be observed 
using two-tone spectroscopy. The transmission through 
one of the resonances inside the band is measured, while 
sweeping a second probe tone that is applied by an addi¬ 
tional microwave source (Agilent 8257D) through a power 
combiner (Mini-Circuits ZFRSC-183-S). This detection 
method is based on the nonlinear interaction between 
the modes. We plot the phase of the transmitted signal 
through the j = 2 resonance as a function of the probe 
frequency in Figl3l We observed the first 14 modes out 
of the 199 predicted modes of the chain. 

In figure IH we show the dispersion relation (frequency 
as a function of the mode index) extracted from the two- 
tone spectroscopy in figure [3l As will be shown in detail 
below, the observed frequencies are in good agreement 
with the theoretical prediction. 
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FIG. 4. Measured frequencies (extracted from the measure¬ 
ment shown in figure [3|) of the resonances of a chain of 200 
SQUIDS compared to the theoretical prediction. The values 
of the ground capacitance as well as the coupling capacitors 
were used as fit parameters in order to find best agreement 
with the data. The highlighted area corresponds to the band 
in which a direct spectroscopy is possible with our measure¬ 
ment setup. 


III. THEORY 


In this section, we will analyze our experiments based 
on the circuit model depicted in figure [5l It consists of a 
chain of N junctions (the SQUIDs are modeled as single 
junctions) with Josephson energy Ej and capacitance C. 
The ground capacitance of the superconducting islands is 
Cq. The chain is coupled to the outside world via input 
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and output capacitors Cin and Cout at the two ends of 
the chain. The corresponding Lagrangian reads 




N 


n=l 

N-1 


1 . 


N-1 


^ ^ ^ (^n+1 ^n) C ^ ^ -S'j COs(0^_|_i 4^n)'i 


n=0 


n=0 


( 2 ) 


where we denoted the superconducting phase of the n- 
th island as 0 ^, and the corresponding node flux = 
h(j)nl2e. In the absence of a voltage bias, the input 
and output capacitors are grounded; we then set = 
^out = O 5 the reference phase for ground being chosen 
equal to zero. Below, we will first analyze Lagrangian ([ 2 ]) 
by linearizing the Josephson term, i.e. expanding it up 
to the second order in (0n+i — Then we will study 

the effect of weak nonlinearity perturbatively by includ¬ 
ing fourth order corrections oc (0n+i — 0n)^- 


Cjn 

^in 








■'out 


$ 


out 


Co 




FIG. 5. Circuit diagram considered for the derivation of 
the Hamiltonian. A series array of Josephson junctions with 
Josephson energy Ej, capacitance C and ground capacitance 
Co coupled to the external circuitry via the capacitors Cm 
and Cout- 


A. Linear modes 

In the linear limit, Lagrangian ([2]) of the system can 
easily be written down as a sum over the islands of 
quadratic terms only 


N-1 N-1 

+ E 2(^"+i - - E 

n=0 n=0 

(3) 

where we introduced the Josephson inductance L = 
{h/2e)‘^/Ej. In the last line, we introduced the matrix 
form with the flux vector ^ and its transpose defined 
as = (<I>o, 4 >i,..., 4 > 7 v) and the capacitance matrix C 
and the inverse inductance matrix L~^. 


From Lagrangian © the Hamiltonian can be derived 
via Legendre-transformation 


Ho =Q^I - £ = ^Q^C-^Q + (4) 

where the components of the charge vector are defined 
by 

The Hamiltonian (|4]) can be quantized and expressed in 
terms of the usual bosonic creation and annihilation op¬ 
erators. Details are presented in appendix [Al The Hamil¬ 
tonian 0 ] then reads 


N 

Hq = Tiuoj 
i=i 


a]aj 


(5) 


where the eigenfrequencies are given by the eigenvalue 
problem 


C-^/^L-^C-^/^Vb- = (6) 

with = C~^. The eigenvectors ipj of the 

matrix are related with the eigenmodes 

of the Hamiltonian (|4]) via 

* = + (7) 

The eigenvalues ujj of the Hamiltonian © are very 
sensitive to the values of the coupling capacitors Cin and 
Cout- This is because these capacitances determine the 
boundary conditions for the allowed eigenfunctions. To 
illustrate this, we plot the eigenfrequencies of the lowest 
modes as a function of the input capacitance in figure [ 6 l 
We keep a fixed ratio between output and input capaci¬ 
tors; the other parameters were chosen according to the 
sample parameters listed in Table HI We see that, upon 
increasing the eigenfrequencies decrease monoton- 
ically. Indeed, when Cin = Cout = O 5 fhe chain is ef¬ 
fectively isolated; as a result the allowed eigenfunctions 
satisfy the zero-current boundary condition <I>o — 4>i = 
4 > 7 v-i — 4 > 7 v = 0. The corresponding N — 1 allowed modes 
are given by 4>n = AkCos[k{n — 1/2)], where the con¬ 
stant Ak fixes the normalization and the dimensionless 
wave vector k = jtt/ {N — 1) with j = 0,1,..., A — 2. In 
the opposite limit of large C^n, and hence Cout^ islands 
0 and N are perfectly coupled to the outside electrodes 
connected to ground, thus the boundary condition reads 
4>o = 4 > 7 v = 0. The corresponding N — 1 modes are 
now given by = Bf^sin{kn), where the constant 
fixes again the normalization and the dimensionless wave 
vector k = jTrjN with j = 1, 2,..., N — 1. The modes fre¬ 
quencies follow from the usual dispersion relation, valid 
for homogeneous chains^ 


1 — cos k 


— ^p\ 


1 — cos k + Co / 2 C 


(8) 
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Their limiting values are indicated in Fig. [6l both for 
Cin = 0 and for Cin = oo. The numerical solution for 
finite Cin smoothly interpolates between these limiting 
values. The gradual change of the nature of eigenfunc¬ 
tions as Cin is increased is shown in Fig.[3 One sees that 
the eigenfunctions indeed have zero slope at the entrance 
and exit of the chain for small whereas they satisfy 
the zero-phase boundary condition for large Cin- 


itance Co and the coupling capacitances Cin and Cout 
in order to find best agreement with the experimentally 
observed resonance frequencies. The result is shown in 
Fig. m where we plot the observed frequencies as a func¬ 
tion of the mode index j. The observed frequencies are 
in good agreement with the theoretical prediction. 


B. Weak nonlinearity 



FIG. 6. Eigenfrequencies as a function of the coupling capac¬ 
itor. We fixed Cout = 3Czn. All other parameters were taken 
from the sample parameters listed in Table [H The eigenval¬ 
ues were calculated using Eq. (|6|). The black vertical line 
indicates the coupling capacitances realized in the experi¬ 
ment. The solid dots indicate the analytical values for small 
{Cin = 0) and and large {Cin oo) coupling capacitors. 

When comparing the theory with the experimental 
data, we have adjusted the values of the ground capac- 


We now include the weak nonlinearity, namely the 
quartic term of the Josephson energy —Ej{(l)n-\-i — 
(/)n)^/24 in the Hamiltonian and treat it as a perturba¬ 
tion to the linear Hamiltonian The perturbative 

approach is justified if the mode frequency ujj <C Ej/Ti 
which is the case for the low frequency modes of our sam¬ 
ple. 

Using the mode expansion ©, the modified Hamiltonian 
Hnl = Hq ETJnl can be found, with a nonlinear contri¬ 
bution Uatl, 



> ~ E 

3 

3 

- E 

— Kj^a^^ajO^^ak 

j,k 

(j 7^ k) 


- E 


j,k 

U + k) 


- E 


j, k, 1 



{j^k^ 1) 


(9) 



EIG. 7. Space dependent part of the eigenfunction for the 
j=2 mode for different values of the coupling capacitance Cin- 
{Cout — 3Gin) 


The coefficient Kjj is called self-Kerr coefficient and the 
corresponding self-Kerr term of the Hamiltonian (j9j) gives 
a frequency shift of the frequency of the mode j that 
scales linearly with the photon number a^-aj in the mode 
j. At higher photon number, this term can give rise 
to photon blockade^liSS and bistability of the resonator. 
The term containing the cross-Kerr coefficient Kjk causes 
a shift of the frequency of the mode j that depends lin¬ 
early on the photon number in the mode k. As will be 
discussed in more detail below, the cross-Kerr effect can 
be used to probe the photon numbers. The cross- and 
self-Kerr coefficients calculated with Eqs. m and m 
for the parameters of the sample studied in this paper 
are listed in Table [III Finally, the two terms containing 
the coefficients Qki give rise to a coupling between two 
different modes that depends on photon numbers in a 
third mode. These terms can also give small corrections 
to the self- and cross-Kerr coefficients in second order. 
These corrections are small by a factor i^/A, where K 
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mode index 

2 

3 

4 

2 

2.8 ±0.4 MHz 

6.6 ±0.9 MHz 

9.4 ± 1.4 MHz 

3 

6.6 ±0.9 MHz 

8.3 ±0.7 MHz 

15.9 ± 2.4 MHz 

4 

9.4 ±1.4 MHz 

15.9 ± 2.4 MHz 

16.7 ± 2.5 MHz 


TABLE IL Kerr coefficients for the modes j=2,3 and 4 pre¬ 
dicted by our theory. The errors were estimated assuming the 
maximum errors for C and Cq. 


is the relevant Kerr coefficient and A the difference be¬ 
tween the mode frequencies involved in the second order 
process. 

The nonlinear coefficients are given by 


2nTT^Ejrijjjj 

( 10 ) 

Ahn'^Ejrjjjkk 

( 11 ) 

ihn^Ejrijjki 

( 12 ) 


and the renormalized mode frequency is given by 


^j—^j (13) 

k 


Here we introduced the dimensionless quantities 


vjkip = Yl 

n L 

m 

m 


(14) 


that take into account the spatial variation of the 
modes.2prn,j is the m-th component of the vector de¬ 
fined by Eq. [6l 

The details of the derivation of Eqs. ®, (d and m can 
be found in appendix B. It is noteworthy that although 
the expressions for the Kerr coefficients (Eqs: llQlllll2p 
explicitly contain the Josephson energy Ej^ the Kerr co- 
efficents actually do not depend on Ej. Indeed the mode 
frequencies ujj also contain a factor a/Ej so that the de¬ 
pendence on the Josephson energy cancels out. 


IV. EXPERIMENTS IN THE NONLINEAR 
REGIME 

A. Photon number calibration with the self-Kerr 
shift 

The self-Kerr shift can be used to calibrate the photon 
number in a resonant mode by measuring the resonance 
frequency as a function of the drive power. The resonance 
frequencies as a function of power are obtained by fitting 
the resonance curves for each excitation power. If only 
one mode of the chain is excited, the Hamiltonian Eq. m 
reduces to 


Hsk — ( 


a]aj. 


(15) 


The mean photon number rij of a mode is related with 
the incident power Pin via the relation^^ 


47j,in p 


(16) 


Here 7j,tot is the total decay rate of the mode related 
with the width of the resonance, ^j^in is the decay rate 
through the input capacitor. Pin is the incident power to 
the sample. Equation [16] can be rewritten in terms of the 
input power at room temperature 


rii 




(17) 


where A is the attenuation of the input line of the cryo¬ 
stat, P is the input power at room temperature and aj 
is a numerical factor that takes into account the uncer¬ 
tainty of the attenuation A, reflections at the input of 
the sample as well as the ratio 7j,in/7i,tot- 

We use the self-Kerr measurements to extract the pa¬ 
rameter aj for each of the modes. In figure [8] we plot the 
experimentally observed self-Kerr shifts together with the 
theoretical prediction. 

The values of aj that were used to adjust the photon 
numbers in the modes are shown in the figures. The non- 
uniform variation of aj with the mode index is related 
with spurious resonances in the sample holder. Indeed 
the transmission background (off resonance) is highest 
for the j = 3 resonance and smallest for the j = 2 reso¬ 
nance. 

In the low power region (the few photon limit) the line 
shape of the resonance is altered very little by the pres¬ 
ence of the non-linearity so that the fitting with the Fano 
formula Eq. m is possible. In the limit of high drive pow¬ 
ers the fitting with the linear model Eq. m becomes less 
accurate. The vertical error bars in FiglH^md FiglO] rep¬ 
resent the error of the fit to the resonance curve. 

The perturbation approach breaks down at higher pho¬ 
ton number as the terms KjjUj or KjkUk are no longer 
small compared to the mode frequency ujj . 

For the high excitation amplitudes the power dependent 
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FIG. 8. Experimental self-Kerr shift compared to theory. We 
used a factor aj to adjust the experimental photon number. 
The three graphs correspond to the three modes of the chain 
j = 2 (top), j = 3 (center) and j = 4 (bottom) that lie within 
the bandwidth of the measurement setup. 


frequency shift is reversed so that an increase of the reso¬ 
nance frequency with drive power is observed. The onset 
of this upwards shift is nevertheless visible in figure [8] at 
high photon numbers. 

The errors on the photon number calibration (aj) dom¬ 
inate over the uncertainty of the theoretical prediction 
caused by the uncertainty of the sample parameters (Ta¬ 
ble: HI). aj depends on the range of drive powers that is 
used for the calibration. We estimated the errors of aj 
by performing the photon number calibration for differ¬ 


ent power ranges. 


B. Measurement of the cross-Kerr effect 

With the photon number calibrated we are now able to 
measure the cross-Kerr shifts. For two-tone driving the 
Hamiltonian Eq. (|9]) now takes the form 

Hck = (nojj - - ^Kjkfik^ a]aj. ( 18 ) 

To experimentally investigate the cross-Kerr shifts, we 
measured the transmission through the mode j with the 
network vector analyzer with fixed excitation power, and 
varied the excitation power of the mode k with an ad¬ 
ditional microwave source. As the excitation power of 
the mode j is kept fixed, nj is constant and the self-Kerr 
effect only gives a global offset. To convert the exci¬ 
tation power of the mode k to the photon number, we 
used the value of the ak that we previously calibrated 
with the self-Kerr shifts. The bare resonance frequen¬ 
cies 00 j — ll2Kjjnj were also extracted from the self-Kerr 
measurement. The measurement of the cross-Kerr coef¬ 
ficients therefore does not contain free parameters. 

The measurements and the predicted shifts are plotted 
in Fig. [9l The vertical error bars correspond to the er¬ 
rors of the resonance fit whereas the horizontal error bars 
represent the estimated errors of the photon calibration. 
We observe frequency shifts that are of the same order as 
the shifts predicted by the theoretical model. The K ^2 
could not be extracted from our measurements. 

The cross-Kerr shifts were measured with a fixed drive 
frequency for the mode k. Although the self-Kerr shift 
of the mode k is smaller than the line width of the reso¬ 
nance, this can cause deviations of the linear dependence 
of the photon number with drive power (as the drive can 
become slightly off-resonant) and thus to deviation from 
linear frequency shift for the mode j. This could ex¬ 
plain the curvature observed in some of the cross-Kerr 
measurements. Errors in the chain parameters C and Cq 
affect self- and cross-Kerr coefficients in the same way. 
As our method of using the self-Kerr shift to calibrate 
the photon number and the measuring the cross-Kerr 
shift probes the ratio between self- and cross-Kerr co¬ 
efficients our result is only weakly affected by such er¬ 
rors. We have also investigated the effect of disorder on 
the mode structure of the chain (not shown)^. The low 
frequency modes turn out to be robust against disorder 
as these long wave length modes average over a large 
number of junctions. For the same reasons also their 
Kerr-coefficients are immune to disorder. 


V. SUMMARY & CONCLUSION 

We have investigated extended plasma resonances in 
chains of Josephson junctions. We have measured a 
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j = 2 


j = 2 



j=4 


mean photon number 

j=4 




FIG. 9. Cross-Kerr shift Kjk between the j — 2^ j — and j — ^ mode as a function of the photon numbers in the mode 
k — 2^ k — and k — A. The photon number in the mode k was calibrated with the self-Kerr shift. For two drives the Kerr 
Hamiltonian takes the form Hck = {hujj — ^Kjjfij — ^KjkUk) 


SQUID chain containing 200 junctions and were able to 
measure the 14 lowest resonant modes of the chain in 
a two-tone scheme. The observed resonance frequencies 
are in excellent agreement with the theoretical prediction 
that we obtained by diagonalizing the Hamiltonian in the 
linear limit. Three of these plasma modes lie within our 
measurement band and can be accessed directly in the 
experiment. 

We have presented a derivation of the Kerr-coefiicients 
of the extended plasma resonances of Josephson junction 
chains. The non-linearity of the Josephson inductance 
was therefore treated in perturbation theory. To compare 
the theoretically predicted frequency shifts to the shifts 


observed in the experiment, we use the self-Kerr effect 
to calibrate the photon numbers of the plasma modes as 
a function of the drive power. Then we compared the 
cross-Kerr shifts with the theoretical prediction and find 
reasonable agreement. 

When studying the superconductor insulator transition 
in Josephson junction chains, the low frequency plasma 
modes constitute the electromagnetic environment of the 
junctions. This environment provides quantum fluctua¬ 
tions that enables quantum phase slips on the junctions. 
It is the interplay between the plasma modes and the 
quantum phase slips which depends on the ratio Ej j 
that characterizes the superconductor-insulator transi- 
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tion. 

Here we have studied for the first time the non-linear 
interaction between the plasma modes. The signature 
of the superconducting-insulating transition on the mi¬ 
crowave transmission spectrum of a JJ chain is currently 
under theoretical investigation. The theoretical descrip¬ 
tion as well as the experiments presented in this work 
focus on uniform Josephson junction arrays but an ex¬ 
tension to non-uniform arrays is straight forward. Indeed 
recent works^S discuss mode engineering of the resonant 
modes of a chain by varying the junction parameters. 
This approach might be extended to also engineer the 
non-linear interaction between the modes and thus tailor 


an active medium for the use in travelling wave paramet¬ 
ric amplifiers^i^. 
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Appendix A: Linear Hamiltonian, eigenvalues and 
eigenfunctions 


In this Appendix we give a rigorous derivation of the 
theoretical model presented above. Let us start from the 
Lagrangian Eq. (|3|) which can be written in the compact 
form 

n = (Al) 

with the flux vector 

= (i>o, $1 ,..., i>Ar), (A2) 


and the capacitance and inverse inductance matrices de¬ 
fined as 


( a -\- Gin 


C = 


-C 

0 


-C 0 
Co + 2C -C 0 
-C Co + 2C -C 

0 


0 




\ 


/ 

(A3) 
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and 


Appendix B: Nonlinear Hamiltonian 


L-i = 


/I 


L 

0 


-1 

L 


-1 

L 


-1 

L 


-1 

L 


(A4) 


The momentum vector associated with the flux vector 
has the units of a charge and is given by 

Q = Ci. (A5) 

With this the Legendre transformation gives the Hamil¬ 
tonian O- 

In order to derive the eigenvalues and eigenfunctions 
of this Hamiltonian, let us write the total energy of the 
system as 


where 

0 = (A7) 

and = C~^. A Fourier decomposition gives 

Ht) = E . (A8) 

3 

2 pj and ujj satisfy the eigenvalue equation 

g-l/2£-lg-l/2^_^2^^._ (A9) 

Inserting eqs. (|A8p and (|A9p back into eq. (|A6p we find 

E = (AlO) 

3 


The energy of the normal modes and thus also the Hamil¬ 
tonian can be expressed in terms of creation and annihi¬ 
lation operators Hq = Tiujj{o}-aj -h 1/2) where 


ipj = 



(All) 


The flux operator for the fluxes on the islands can be 
obtained by inverting Eq. (113 and inserting Eqs. dm, 
(ESI) and (lAnTi 


We include the lowest order nonlinear term of the 
Josephson energy, -Ej - 4>nt as a per- 

turbation to the linear Hamiltonian Eq. &■ We expand 
this nonlinear term using the eigenfunctions Eq. (jA12p . 
As a result we obtain the weak nonlinear contribution 


1 / 9 \ N—1 

^ A n=0 


1 IbTT^ 


N-1 

24 ^ 

0 n=0 




m,j 








3 ^ 




j,k 

O' ¥= 'm) 


2h TT Ej \ Vjjjk 

3$oC2 ^ LOj^/uJjLOk 

j,k 

O' ^ k) 

h^TT^^Ej 


(®i + ®i) + “I) 


$4^2 


E 

j, k, I 

0 ^ M 0 


Vjjkl 


Uljy/UlkWl 


(dj dj) 


h^TT^^Ej 

6$^C2 


X ^dfe -h d[) (^d; + d^) 

^ ^CUjOJkOJiOJm ^ •''' 

j,k,l,m 

U ¥=k^ly^m) 

X ^dfc + d^) ^d; -|- d|) (dm + d|)j) , 

(Bl) 


where we introduced the superconducting flux quantum 
$0 = h/2e. Equation (|9]) in the main text is then ob¬ 
tained by using 


(^dj -h a] j = baja^djdj -h baja^ + 3, (B2) 


$ = Ea-i/Vi 

3 



(A12) 


(% + + ®l) 

= 4d]djd^dfc -I- 2d]dj -t 2d[dfe -|- 1, (B3) 
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(^dj + a]) (dk + 

= {3d]dj){d]dk + 4«i) + 34«i> (B4) 

(«i + Oj) («fe + «l) 

= (2a]aj + l)(afeaz + djdk), (B5) 


and 


(^dj + a]) (^dk + 


(^di + dl'j {dm + dL) = ajdkdjdr, 


+ . 


(B6) 


